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Abstract 



We report our finding of a two scale cut off structure in the in- 
frared nonperturbative dynamics of Pade couplant QCD, in the flavor 
states Nf < 8. We argue that these two NPQCD momentum scales 
Q c (Yi) and Qc^Xt) can be identified as the scales of onset of chiral 
c~| , symmetry breaking and quark confinement in QCD, and as the cut 

off boundaries of a strongly interacting infrared quark gluon phase of 
Quantum Chromodynamics, intermediate between the hadronic phase 
of QCD and the weakly coupled perturbative QCD regime. We inves- 
tigsted the pattern of variation with flavor number Nf < 8, of these 
domain boundaries of the intermediate nonperturbative infrared QCD, 
and found that the two scales are correlated, with a correlation factor 
QcOfy/QcOfl) that rises to a peak at Nf = 2 from Nf = 0, but falls 
off rapidly to zero for 4 < Nf < 8. We concluded firmly that dy- 
namical chiral symmetry breaking and quark confinement while being 
two distinct QCD phenomena caused by two independent component 
QCD forces Y\ and Y2, are nevertheless closely related phenomena of 
infrared nonperturbative QCD dynamics. Their correlation leads us 
to a finding that quark confinement is most favored in the Nf = 2 
flavor state of QCD, but becomes rapidly less probable for Nf > 4, 
this finding being exactly as one observes in nature. 
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1 INTRODUCTION 

Quantum Chromodynamics (QCD) is known to have three main regimes 
or phases: (1) the perturbative QCD (PQCD) regime at high momentum 
transfers Q where the color force dynamics is weak and tends to zero for Q — ► 
oo; (2) the hadron regime at very low momentum transfers (the far infrared) 
where the color force has already grown so strong that the quarks and gluons 
are confined inside hadrons and cease to be the observable degrees of freedom 
of QCD dynamics. (3) Then there is the intermediate (infrared) momentum 
transfer regime we shall here denote by CQCD, where the quarks and gluons 
are still the primary physical degrees of freedom of the QCD dynamics, but 
have their coupling constant sufficiently large that all the QCD processes in 
this regime are nonperturbative. These nonperturbative processes include 
in particular chiral symmetry breaking and the confinement of quarks and 
gluons into hadrons. As conceived, the CQCD regime of QCD necessarily 
interfaces with the regime of applicability of Chiral Perturbative theories 
and Effective Lagrangians Jl|, ||, Q of QCD, where Nambu-Goldstone bosons 
(pions, kaons..) formed from nonperturbative chiral symmetry breaking, are 
explicitly highlighted as the effective degrees of freedom in terms of which to 
formulate effective Lagrangians of QCD dynamics in this low energy region. 

The exact momentum scale points and critical couplings at which the 
CQCD regime sets in from the PQCD end, and later cuts off at the hadron 
end, are not known or measured, although a number of phenomenological 
studies particularly of chiral symmetry breaking, quark antiquark conden- 
sates, and quark confinement ||, ||, || |7J, have placed some bench marks on 
the onset and cessation of this CQCD regime of Quantum Chromodynam- 
ics. Additionally, recent heavy ion experiments at CERN and Brookhaven 
(RHIC), approaching the same CQCD regime of QCD from the hadron end 
(at high temperatures), have confirmed the qualitative existence of this third 
phase of QCD called there Quark Gluon Plasma (QGP) phase, and begun to 
yield some quantitative information on the points of onset of deconfinement 
and chiral symmetry restoration in high temperature QCD. Finally, Lattice 
Quantum field theory simulations of QCD [§, || have several hints of chiral 
symmetry breaking and quark confinement phases in QCD, and have placed 
their own bench marks on the domain of the QGP intermediate phase of 
QCD. 

These known features of QCD, in particular the three phase structure 
of QCD, lead us to consider seriously the recent finding in our paper |TTJ, 



based on Pade improved PQCD, that the same low energy region of QCD is 
governed by four component coupling constant solutions 01,02,03,04, (with 
Oj = a^/V), in which a\ and a 2 account for any pole structure or else coupling 



constant freezing in the low energy region of QCD, while 03 and a 4 provide a 
two scale cut off structure of PQCD from the NPQCD infrared region. The 
possibility then arises that these two IR cut off scales a 3 and a 4 observed in 
our Pade couplant QCD, may be identified as the onset points or boundaries 
of the CQCD regime of QCD, and in particular as the critical points of chiral 
symmetry breaking and quark confinement in QCD. This possibility is what 
we examine in this paper, matching our observed two scale Pade infrared 
NPQCD coupling constant values and features (a 3 , a 4 ), with various existing 
phenomenological bench marks and features of low energy QCD dynamics. 
Because no such two scale infrared cut off structure has previously been 
observed or reported in perturbative QCD studies with renormalization group 
equations, in the infrared region, our work appears to us an interesting new 
finding. 

We present our study and results as follows. In section 2, we recall some 
essential equations and features of our optimized [1|1] Pade QCD from our 
earlier paper |10|| , noting in particular the observed two scale Q c (Yi) and 



QcOty infrared cut off structure in the flavor states Nf < 8 of our opti- 
mized Pade couplant QCD. Because this earlier observation was based on 
the one physical observable or effective charge R e + e -(Q) and its renormal- 
ization scheme invariant P2, we devote sections 3 and 4 of the paper, to con- 
siderations of what these Pade couplant features look like, when other QCD 
physical observables, both timelike and spacelike, are used besides R e+e -(Q). 
Satisfied that the two scale infrared cut off feature persists in all QCD phys- 
ical observables, timelike and spacelike, we devote section 5 of the paper, to 
comparing this two scale feature with the known phenomenological results 
and bench marks of chiral symmetry breaking and confinement in QCD. 

Prompted by the good agreement we found, we devote section 6 of the 
paper to a consideration of the flavor dependence of our infrared two scale 
structure and the implication of this for quark confinement in QCD. Section 
7 considers similar implication for the question of infrared fixed point in low 
flavor states of QCD. The paper concludes with a brief summary in section 
8. 



FEATURES OF THE OPTIMIZED [1|1] PADE 

QCD WITH R e+e -(Q) OBSERVABLE. 



In our earlier paper [ IQfl , we showed that the optimized [1|1] Pade improved 
PQCD couplant equation bifurcates into two independent Y\ and Y 2 couplant 
equations given by 
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and 



D = 8c 2 - 4p 2 . (6) 

Here a = a s /n, is the optimized couplant solution of eqn. ([[]) or (D), 
while a s is the corresponding QCD coupling constant. We shall from here 
onwards drop the bar over both quantities purely as a matter of notational 
convenience but it is understood that all our computations and plots below 
refer to a or a s /n. The quantities b, c, C2(MS).. are perturbative QCD beta 
function coefficients having their usual values O, [L2|, [13], [IJJ: 
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c 3 (MS) = fa/fa, where 
& = 114.23033 - 27.133944iV) + 1.5823791./V 2 + 5.85669582 x 10- 3 N'f (10) 



In the optimization process [10] , the coefficient c 2 as a variable satisfies a 
quadratic equation given by: 



c 2 - (16c 2 )c 2 + (4c 2 p 2 + c 4 ) 
whose solutions are: 
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C2 ^c 2 = 3c 2 ±c^D = 2c (- c ±-Vd) (12) 

It follows from this that if the quantity D is negative, c 2 is not real and 
our Y\ and Y 2 equations have no real solutions unless we work with \D\. Since 
D = 8c 2 — 4p 2 , this viability of of eqn. (|TT|) depends in turn on the sign and 
magnitude of the renormalization scheme (RS) invariant p 2 . 

These RS invariants pi, p 2 , P3-... arise in general JTJ, [H], [TJ], [Tj| from a con- 
sistency requirement between the PMS/ECH constraints on a perturbatively 
truncated generic physical observable 

R {n \Q) = a(l + ri a + r 2 a 2 + ...) (13) 

and the basic renormalization group equation (RGE) for the couplant a. 
Explicitly, the RS invariants are given by: 

Pl (Q,N f ) = T-n (14) 

P2(N f ) = r 2 + c 2 - (n + -c) 2 (15) 

or simply 

P2(N f ) = r 2 + c 2 -r 2 -cri (16) 

Their explicit numerical values and sign then depend on the particular 
physical observable and its loop order coefficients (ri, r 2 , T3....) we use in eval- 
uating the RS invariants. In the particular case of the R e+e -(Q) observable 
used in paper [jID| , the p 2 values are large and negative and eqns ([!]) and (||) 



have four crossing point solutions: Y\ = oi, a 2 , 03 and y 2 = 04. 



Our finding in |10] is that solutions a\ and a 2 determine an infrared 
attractor pole behavior of Pade QCD for Nf < 8, while the same ai,o 2 
solutions determine an infrared fixed point frozen couplant behavior for Nf > 
9. On the other hand, the 03 and 04 solutions determine the points of cut off 
of PQCD from the infrared region. In the particular case of Nf < 8 where 
03 and a 4 remain perturbative long after solutions a\ and a 2 have turned 
away from the infrared region, shown typically in fig. 0, the ultimate cut off 
points Q c (Yi) of 03 and Q C (Y 2 ) of 04 can become phyically meaningful, being 
two possible scales of NPQCD infrared dynamics that our Pade couplant 



QCD predicts. In this respect, we found earlier in paper |TI| that while 03 
cuts off typically at a critical momentum Q c (Yi) ~ 1.50 GeV, the a^ cuts 
off typically at Q C (Y 2 ) w A QCD « 300 MeV, shown in figs, g and |. The 
suggestion is that these two cut off points (and their corresponding critical 
coupling constants) may be marking the domain of CQCD, such that we can 



write: Q C (Y 2 ) < Qcqcd < Qc(Yi), with Q C (Y{) = A x being point of onset of 
chiral symmetry breaking, while Q C (Y 2 ) = A con is the critical scale point of 
confinement and phase transition into hadrons. 

Before however drawing this conclusion, we need to examine two points. 
The first is to determine that the above 03,04 optimized Pade two scale 
infrared couplant feature persists whatever the physical observable and its 
renormalization scheme invariants pi we use, particularly spacelike and time- 
like physical observables. The second point is to determine that our numerical 
values of Q c (Yi) and Q C (Y 2 ) from various physical observables and p 2 values, 
timelike and spacelike, agree with and are consistent with what one has come 
to expect from various existing phenomenological studies of chiral symmetry 
breaking and confinement in QCD. We now consider these two issues. 

3 VARIOUS SPACELIKE AND TIMELIKE 
QCD OBSERVABLES AND THEIR IN- 
VARIANTS p h p 2 . 

As already indicated above, while optimization and the entire effective charge 
formalism [|15| , |19| , |20| , [2l|j achieves the result of renormalization scheme (RS) 
independent perturbative QCD, the prize paid is that a QCD (or Pade QCD) 
coupling constant becomes defined and tied to a particular physical observ- 
able and its renormalization scheme invariants pi. If therefore we observe a 
particular feature of such optimized QCD (or optimized Pade QCD) with re- 
spect to one physical observable and its RS invariants, we cannot be sure that 
this feature is universal and holds for other physical observables and their 
respective RS invariants. This substitute ambiguity applies to any finding 
in these optimized QCDs (or optimized Pade QCDs), of couplant freezing 
as well as to the two scale a 3 , a 4 infrared cut off structure discussed above. 
The closest we can go towards ameliorating the ambiguity is to compute 
an observed couplant feature using several different QCD observables, both 
timelike and spacelike, and to accept as reasonably valid and universal, only 
those couplant features that appear reasonably independent of any particular 
physical observable or effective charge used in computing them. To achieve 
this, we first compute here the p 2 values of a number of spacelike and timelike 
QCD observables, whose perturbative coefficients r\,r 2 ,r^.... are known, at 
least up to NNLO loop order. We then test out their optimized [1|1] Pade 
couplant structures in the infrared, and compare with our earlier findings in 
paper JTUJ based on R e+e _(Q) observable. 



The following are the spacelike and timelike QCD observables (effective 
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Figure 1: A plot from the R e+e _(Q) case, showing that Pade couplant solu- 
tions 03 and 04 in general remain perturbative in the region Q < Qmin and 
for all flavors Nf < 8, except close to their Landau type cut off point. 
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Figure 2: The infrared cut off points of Pade couplant solutions a^ and a^ for 
various flavors Nf < 8, computed with respect to the R e+e _(Q) observable 
(effective charge) 
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Figure 3: Further view of the infrared cut off points of Pade couplant 
solutions CL3 and a 4 for various flavors Nf < 8, computed with respect to the 
R e+e ^(Q) observable (effective charge) 



charges) we consider: 

1. Bjorken sum rule for deep inelastic scattering (DIS) of polar- 



ized electrons on polarized nucleons [p2 



l dx{gT{x,Q 2 )-gT{x,Q 2 )} = \ 
o 3 



9a 



9v 



,(e) 



1-RfAa)] (17) 



?(e) 



where R B j(a) is a QCD spacelike observable given perturbatively by: 



Rgj(a) — a(l + ki<2 + k 2 a 2 + 
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Direct loop computations in the minimal subtraction renormalization 



scheme MS p3[ give 

fci = Nt 

12 3 s 



k-, 



13841 44 , 
216 9 s 



55 . 10339 61 5 

yCs - ^/(^age- + 54 C3 " 3 Cs) 



115 ,t2 

648 f 



(19) 



Taking (3 = 1.2020569031; Cs = 1.0369277551, and combining with 
eqns. (0), (||), @ and (0) we compute the values of p 2 shown in 
Table [J 



2. Bjorken sum rule for DIS of neutrinos on nucleons [22] 



dx{F?{x,Q 2 ) - Fr P (x,Q 2 )} = 1 - ^(41) 
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where R^j(a) is a QCD spacelike observable given perturbatively by: 
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J Rg(a) = a(l + / 1 a + / 2 a 2 + 



(21) 



Direct loop computations in M5* [0] give: 

h 
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Again one computes the corresponding p 2 values shown in Table [I], 
from eqn. Qlfip. 
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3. Gross - Llewellyn Smith (GLS) sum rule |25 



dx{F?{x, Q 2 ) + F?{x, Q 2 )} = 6[1 - R^ s 



M 



(23) 
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where RGLs( a ) * s a QCD spacelike observable given perturbatively by: 



R GLs( a ) 



o(l + SiO + S2CI + 



(24) 



Direct loop computations in MS [23, 26] give: 



55 1», 
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10009 91 5 AX 115. 
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(25) 



The corresponding p 2 values we computed from eqn. (JI 

shown in Table 0. 

4. Total hadronic cross section R(s) in e + e annihilations 



are 



i?(s) = 3^gj(l + i? e+e -(a)) 
/ 



(26) 



where R e + e - (a) is a QCD timelike observable given perturbatively by: 



Re 



all + rio + r2<i + 



(27) 



with coefficient values given in MS given by [271, pB] 



31 



ri = 1.9857 -0.1153JV/ 

T2 = — u.uouo — ±.Z,UU±iVy 



«.Ui368 ] .200] Xf - 0.0052iV? - 1.2395 „„£; (28) 
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3EQ 2 



/ 



Ae+e-) 



The corresponding p 2 = P2 values were already computed from 

eqn. (|15| ) and used in our paper |10[. They are reproduced here in 
Table [l] for direct comparison. Formula eqn. (|15|) tends to make the p 2 
values larger and more negative at the higher flavors Nf, than eqn. (|i~6p. 
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5. r Lepton hadronic decay ratio R T [p2l B3, [3 



„ T(t — *■ z/ T + hadrons) , . 

^ = l r , — — v^ 29 

This has a timelike QCD component observable R^ ert given perturba- 

tively by 
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(30) 



where 



Ft = (2N f - 33)/6 = -b 

F 2 = (19N f - 153)/12 

F 3 = 1.9857- 0.1 153^/ 

F 4 = -6.6368 - 1.2001iV> - 0.0052iVj (31) 

and V u d and V us are usual flavor mixing (CKM) parameters we do not 
however need here. The p 2 values computed from eqn. (|T6|) are shown 
in Table [3]. 



6. Higgs hadronic Decay Width p5|, |36 



T tot (H -> hadrons) = ^^ £< R s ( s ) (32) 

V 2 ^ Nf f 

where R s (s) is a QCD timelike observable given perturbatively up to 
NNLO by p 



R s (s) = 3 {l + 5.66667a +(35.93996 - l.35865N f )a 2 

+ (164.13921- 25.771197V/ H- 0.258974iVj)a 3 } (33) 

or writing 

R s (s) = 3(1 + R(s, a)) with R(s, a) = o(l + ha + h 2 a 2 + ) 
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we have : 

35.93996 - 1.35865^/ 



h 



5.66667 



(164.13921 - 25.77119iV f + 0.258974iV?) , N 

h 2 = - -i f - (34) 

5.66667 

The corresponding p 2 = p 2 values computed from eqn. ( |l~6|) are shown 
in Table § 



4 PADE COUPLANT FEATURES FOUND 
WITH SPACELIKE AND TIMELIKE OB- 
SERVABLES 

Given now the above p 2 values for different observables shown in Table [I], 



we use the same computational procedures described in paper [|10| , to as- 
certain the multicomponent couplant structures of eqns. ([[]) and (0) with 
respect to each of the above observables and their p 2 values. Our findings 
and comparative features are as follows: 

1. First we note from Table |l] that the spacelike observables Bj(e), Bj(z/), 
GLS(z/) have positive p 2 values for Nf < 4, except for a small negative 
p 2 = —0.3223 at Nf = 4 in the GLS case. Correspondingly, we find that 
the quantity D of eqn. (|) is negative for Nf = 0, 1, 2 in these spacelike 
observables. In these negative D cases, eqn. ( |Tl~D has no real c 2 solution, 
and progress is possible only if we use J \D\ in eqns. (^) and (^). For 
Nf = 3,4, the quantity D is positive even though p 2 remains positive. 
Subject to these observations, equations (f|) and (|5|) apply exactly as in 
the timelike cases, and our observed couplant features should become 
comparable in all spacelike and timelike cases. 

2. Computing the crossing point features of eqn. (P for Nf > 3, we find 
exactly the same triple point Y\ = ai,a 2 ,a^ solution for the spacelike 
observables as for the timelike observables. This is shown in fig. H 



and agrees with the features we found in paper [|10[ . The same triple 
crossing point solutions Y\ = ai,a 2 ,ct3 are also found for Nf = 0, 1,2 
flavor states of the spacelike observables, provided we use |D| in eq. (U). 
They all give similar Y\ spiral structures with parameter values shown 
in Tables 0, 0, and £|, to be compared with similar structural parameters 



given in paper [TQ] for the R e+e _(Q) observable. 
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When we compute next the crossing point features of Y 2 from eqn. (0), 
we find that in general for Nf < 4 the spacelike observables behave 
differently from the timelike observables. The Y 2 plots of these spacelike 
observables for Nf < 3, show a triple point crossing feature: Y 2 = 
bi,b 2 ,b 3 , much like their Y x = ai,a 2 ,a 3 counterpart. The Nf = 3 case 
of this Y 2 triple crossing feature can be seen in fig. [5| In contrast, 
for Nf > 5, all physical observables, spacelike and timelike, show only 
one Y 2 crossing point solution a 4 seen in fig. |5] and comparable to our 
earlier findings in paper flO |. The above triple Y 2 crossing behavior for 



Nf < 3 in the spacelike cases, will be found later to lead to a peculiar 
infrared cut off behavior of Y 2 relative to Y\ in these low flavor states 
of spacelike observables. 

4. We look closely next at the infrared cut off region of these spacelike 
and timelike observables. Taking first the cases of Nf = 5,6... for the 
spacelike observables and looking in their infrared cut off region shown 
in figs. H and [7|, we see exactly the same two scale cut off structure, 
as that shown in figs. ^ and |3] for the timelike R e+e _(Q) case. The 
structures separate QCD dynamics in momentum transfer (or energy) 
space into three distinct regimes we have labelled in the figures as 
PQCD, CQCD, and hadron QCD. 

5. Exactly this same two distinct scale cut off structure is exhibited for 
all flavor states < Nf < 8 by the r hadronic decay rate and the Higgs 
hadronic decay width observables shown in part in figs. || and |9|. The 
figures compare well with figs. ^| and |3] of the i? e+e _(Q) observable. 

6. The peculiar cases of Nf < 4, for spacelike observables are shown in 
figs. [10] to |1^, to be compared with the same flavor state plots of timelike 
observables shown variously in figs. @, EE H an< ^ i- We see that arising 
directly from the positive p 2 values in these Nf < 4, flavor states of 
spacelike observables, their Y 2 and Y\ curves tend to bunch together at 
the infrared cut off, with the Y 2 appearing to cut off earlier than the Y\ 
curve for Nf < 3, though not for Nf = 4. The bunching together tends 
to reduce the CQCD momentum gap or domain, compared to figs. § 
and [| or figs. |6] to |9]. However, the two scale infrared cut off structure 
still remains a finding, even in these positive p 2 spacelike cases. 

7. Having seen that the two scale nonperturbative infrared structure per- 
sists in all flavor states Nf < 8, and for all observables spacelike and 
timelike, we now consider the paired numerical values of these scales, 
Q C {Y 2 ) and Q c (Yi), and how these values vary from flavor to flavor and 
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from one QCD observable (or perturbative series) to another. There 
are two obvious points at which we can evaluate and compare the 
two scales. We can evaluate the two scales at the point where each 
curve Y 2 and Y 1 independently begins to depart appreciably from its 
smoothly rising asymptotically free PQCD curve, to enter a nonper- 
turbative QCD phase characterised by sudden rapidly rising coupling 
constant a^ and a 3 respectively, in the infrared region. We shall denote 
the values of Q C (Y 2 ) and Q C {Y\) so evaluated at this onset point of non- 
perturbative infrared QCD dynamics, by Q on (Y 2 ) an d Qan{Y\) respec- 
tively The second obvious point where we can evaluate and compare 
the two scales is where each curve independently cuts off at the far 
infrared end, signalled by its coupling constant rising to infinite values 
(a > 30, or a s > 100). We shall denote the two scales evaluted at this 
terminal infrared cut off point for each curve, by Q C (Y 2 ) = QoffOty an d 
Qc(Yi) = Qoff(Yi) respectively. A third obvious evaluation point exists 
in the special cases of Nf < 4 for spacelike observables, where as shown 
in figs. [R] to 12, the Y 2 curve exhibits a sudden turning away from the 



infrared region within the nonperturbative QCD regime, before finally 
rising up to infinity. This Y 2 infrared attractor point, provides a dis- 
tinct coupling constant hierachy point at which we can also evaluate 
and compare the two infrared nonperturbative Y 2 and Y\ scales. 

Regarding the onset point of nonperturbative dynamics, we find from 
our various plots shown in figs. ||, ||, and |] to [l^, that the two curves Y 2 
and Y\ commence their transition into the nonperturbative regime at 
about the same critical coupling constant value, seen from the plots to 
lie in the range 0.18 < a con « a x < 0.28, or 0.56 < a c ° n « a* < 0.87, 
where a con denotes this critical coupling as evaluated on the Y 2 curve, 
while a x denotes the same nonperturbative onset critical coupling as 
evaluated on the Y\ curves. If we look at the Y 2 plots as the more 
sharply focused indicator of this point of coupling constant criticality, 
we find from figs. |3|, § and |9|, that the Y 2 curves for different flavors tend 
actually to pass through a common point that lies around a con ~ 0.25, 
beyond which these various flavor lines cross paths and rise sharply to 
infinity. Basing on this, we can take as the mean onset critical coupling 
constant revealed by our plots for both Y 2 and Yi, to be the value 
a con ~ a x = 0.25, Accordingly, we can evaluate and compare Q on (Y 2 ) 
and QoniYi) at a con w a x = 0.25. The values of the two nonperturbative 
infrared scales evaluated at this point (I) and at the other two strategic 
points (III, and II) stated above, are shown in Tables [^ and || 
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Ignoring a few positive p 2 cases where specifically the Y 2 curve cut off 



widely earlier than the Y\ curve as shown in plots (b) and (c) of fig [LI 
we find one running feature of the paired scales shown in Tables |5] and|5|. 
This feature is that of the two nonperturbative infrared cut off scales 
Qoff{Y 2 ) and QoffiYi), the lower one defined by the Y 2 solution occurs 
at a mean critical momentum cut off value of Q ff{Y 2 ) ~ 296 MeV, well 
within a value we would associate with Aqcd, the fundamental scale or 
constant of QCD. On the other hand, the upper cut off scale defined by 
QoffOi) occurs at a mean value: QoffiXi) = 1-44 GeV. Both results 
are consistent with our earlier findings in paper |TD[ based on timelike 
R e+e -(Q) observable. 

10. A second feature we observe readily from Tables |5] and |5], is that irre- 
spective of the point of evaluation and the particular QCD observable 
considered, the two scales Q C {Y 2 ) and Q c (Yi) are functions of flavor 
number Nf, both changing in value as the flavor number runs from 
Nf = to Nf = 8. Here the spacelike and timelike QCD observables 
exhibit a distinct difference in their Nf < 4 states. Explicitly we find 
that for the timelike observables and for all < Nf < 8, the two in- 
frared scales move progressively apart as flavor increases. This can be 
seen directly in figs. ||, [| || and ||. The same feature holds for spacelike 
observables but only for Nf > 4 shown in figs. ^| and [7]. For Nf < 4, 
Tables |5] and |6|, and figs. [TO] to [1^ show that the spacelike observables 



exhibit a different pattern of flavor variation. With a view to later 
investigating what these patterns of flavor variation of the two nonper- 
turbative infrared scales can mean, we have computed and shown in 
Tables |5] and ^ the ratio quantities Q c (Y 2 )/Q c (Yi) of the two scales 
evaluated at one or the other of the three strategic couplant points 
mentioned above for various flavors and QCD observables. The ratios 
have been denoted I, II, II, depending on the point of evaluation of the 
Y 2 ,Y\ scales. Since this ratio gives a measure of the moving apart or 
the closing in of these two scales, we shall call the ratio a correlation 
factor between the two solutions Y 2 and Y x and the dynamical processes 
they cause or represent. 

11. Summarizing, we can affirm that our curves not only define a three 
phase Pade QCD regime labelled PQCD, CQCD, and hQCD in our 
figures, but give us definite critical coupling constant values and crit- 
ical momenta at which one transits from one Pade QCD regime into 
another. The extent these features and findings from Pade QCD lead 
us to conclude about physical QCD are what we examine next. 
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5 PHENOMENOLOGICAL FEATURES OF 

DYNAMICAL CHIRAL SYMMETRY BREAK- 
ING AND CONFINEMENT IN QCD 

To see that the above Pade QCD infrared critical momenta Q C {Y\) and Q C (Y2) 
are amenable to interpretation as boundaries of the CQCD domain of physical 
QCD, we consider now some phenomeno logical features and attributes of 
dynamical chiral symmetry breaking (DxSB) and quark confinement in QCD 
gleaned from a number of different sources and models, as follows: 

1. Miransky || reviewing a broad spectrum of phenomenogical models of 
low energy QCD, in particular Bethe Salpeter wave function models of 
dynamical chiral symmetry breaking (DxSB) and quark confinement, 
concluded that the momentum region of QCD where DxSB takes place 
can be parameterized phenomeno logically as: 5 2 < Q 2 < A 2 where Q is 
the QCD operating momentum, while 5 and A x are phenomenological 
infrared (IR) and ultraviolet (UV) cut-off momenta respectively, the 
values of which were phenomenologically found to be: 6 > Aqcd ~ 300 
MeV, and A x ps 1 GeV. The delta scale 5 ~ Aq CD was identified with 
the confinement scale of QCD, and was thus found phenomenologically 
to be at least a few times smaller than the DxSB scale. These values 
are to be compared with our findings in Tables |5| and |6| where writing 
Qoff(Y 2 ) < Qcqcd < Qoff(Yi) as A con < Qcqcd < A x , we found 
A con = 296 MeV (mean value), while A x = 1.44 GeV (mean value), 
both for all < Nf < 8 and for various physical observables, timelike 
and spacelike. The agreement is good. 

2. Concerning the QCD critical coupling constant at the onset of DxSB, 



Refs. ||, |37|, |38| found from their phenomenological estimates and ap- 
proximations, the value: a} > ir/A meaning a x > 0.25. The same value 
was arrived at by Higashijima || who found in his effective potential 
approaches, that chiral symmetry breaks down when his parameter 
\(t e ) = 3C 2 (R)g 2 /(4n 2 ) > 1, where C 2 (R) = 4/3 for three color QCD, 
implying a} = g 2 /(A7i) > 7r/4, or a x > 0.25. These phenomenological 
QCD critical couplings for DxSB are to be compared with our exact 
finding in fig. |9| and other plots, that the various flavor lines converge 
and rise nonperturbatively at a con ~ ct x ~ 0.25, both for timelike and 
spacelike observables. The agreement is excellent. 



3. Roberts and McKellar |39| solving numerically the Schwinger- Dyson 



(SD) integral equation for quark self energy (condensate), under vary- 
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ing kernel approximations, obtained the specific values: a} = 0.782 
in one kernel approximaion, and a} = 0.890 in a more exact kernel 
evaluation. Their results then place a x in the range 0.253 < a x < 0.29 
for Nf = 0,4 considered by them, to be compared again with our 
above Pade QCD finding in the general range: 0.18 < a x < 0.28 for 
< Nf < 6. The agreement is again excellent. Roberts and McKellar 
additionally found from their work that D\SB and confinement do 
not necessarily occur together; rather the latter is held to occur sub- 
sequently, close to where the QCD coupling constant grows infinitely 
large as momentum transfer goes to zero. Our figs. [I] to [| and figs. ^ 
to |l^ manifest exactly this two scale feature and a sharply rising cou- 
pling constant to infinity, compatible with infrared slavery confinement 
at very low Q values. 



Earlier work by Atkinson and Johnson |40| established the interesting 



point that there exists a critical value of QCD coupling constant, cor- 
responding to the onset of chiral symmetry breaking, provided that (a) 
there is an infrared cutoff, realized phenomenologically as an effective 
gluon mass, and (b) there is an ultraviolet cutoff linked with a running 
coupling constant (a chiral symmetric phase of QCD). Within their 
own approximations to a fermion propagator Schwinger-Dyson equa- 
tion, they find the critical QCD coupling constant for D\SB to have 
the value a} = 0.91. These findings from the Schwinger-Dyson ladder 
approximant equation agree again quite well with our Pade QCD two 
scale infrared NPQCD structure, and with the fact that our plots show 
a * ~ 0.87 at the upper limit. 



5. In more recent studies, Guo and Huang |4l| using two different ap- 
proaches: (a) a semi-phenomenological self consistent equation to de- 
termine the quark condensate (the order parameter of D\SB) in the 
chiral limit, and (b) the usual Schwinger-Dyson equation (SD) im- 
proved by addition of gluon condensate kernel, obtained a critical QCD 
coupling a} above which chiral symmetry breaks down. Their value is: 
a} = 0.27T, which means a x = a^/ir = 0.2. This again agrees very well 
with our Pade QCD findings seen in fig. [I] and figs. [7| to [10, and our 
specific mean value: a con ^a a x = 0.25. 

6. Approaching our CQCD regime of QCD from the hadron end, in heavy 
ion Quark Gluon Plasma (QGP) studies, Lattice gauge theories and 
simulations || ^ fi3|| , have evidence of the same two scale NPQCD 
structure found from our Pade QCD, and in the correct sequence. They 
find that the chiral symmetry restoration temperature T x is somewhat 
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larger than the deconfinement temperature Td eC on- Put differently, Lat- 
tice QCD studies find that D\SB is characterized by smaller distances 
compared to confinement distances. In our Pade QCD results played 
backwards, we expect our lower cut off momentum Q ff(Y 2 ) = A con 
to identify closely with Td econ as two parameters that mark confine- 
ment /deconfinement in QCD. Correspondingly, we expect our upper 
cut off momentum: Q ff(Yi) = A x to also identify closely with T x as 
two parameters that mark chiral symmetry breaking and restoration 
in QCD dynamics. Our Pade QCD finding that in general A x > A con 
then receives good support from the above independent Lattice QCD 
studies. 

7. Explicitly, Shuryak |7J] argued from a variety of phenomenological ob- 
servations, including power corrections to deep inelastic scattering and 
heavy ion quark gluon plasma processes, that there must exist in the 
infrared region of QCD, a second scale other than the confinement 
scale, and that this second scale must be connected with chiral sym- 
metry breaking. He estimated the hierachy of the two scales to be: 
A x sb ^> A con ~ Aq CD , and that their ratio is of the order 3-5, mean- 
ing reciprocal ratios 0.33 - 0.20. Such two scale structure and hierachy 
in the infrared region of QCD is what we find here from our Pade per- 
turbative QCD analysis shown variously in figs. [I] to [5] and figs. ^ to [T^. 
Our Qc(Y2)/Q c (Yi) ratios shown in Tables | and || for our two infrared 
scales are well within the estimates of Shuryak, depending on flavor. 

8. In the context of effective Lagrangian theories of the strongly interact- 
ing regime of QCD where chiral symmetry breaking and quark confine- 
ment are believed to occur, Nambu and Jona-Lasinio (NJL) || were 
the first to introduce a viable effective Lagrangian model of hypothet- 
ical four fermion interactions, in which chiral symmetry breaking with 
pions as Goldstone bosons, was fully realized. The QCD scale at which 
the model worked entered the NJL effective Lagrangian theory as an 
upper cut-off momentum and was explicitly found to be A xSB ~ 1 
GeV. Additionally, Nambu and Jona-Lasinio found that their effec- 
tive theory works in a gap regime of QCD lying between the scale 
A x sb ~ 1 GeV and some lower hadronic scale (~ 300MeV) believed 
to be the onset point of quark confinement into hadrons. What we 
now observe is that this NJL model chiral symmetry breaking scale 
A x sb ~ 1 GeV, as well as the gap regime stretching to lower confine- 
ment scale, are exactly comparable with our finding from our Pade plots 
that, 296MeV < Qcqcd < lAAGeV , the implication of all this being 
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that we identify our Pade CQCD domain with the NJL gap regime of 
Goldstone physics. 

9. Another striking phenomenological model buttressing our Pade case is 
that of Manohar and Georgi || who not only found the phenomenolog- 
ical necessity for the two scale (A x sb, A con ) Goldstone physics structure 
and an intermediate infrared QCD regime bounded by A con < Q < 
A x sb with A x sb 3> A con « Aqcd, but actually proposed the use of the 
ratio of the two scales as a natural expansion parameter in chiral effec- 
tive Lagrangian theories of low energy QCD. Their two scale structure 
has numerical values A x sb ~ 1 GeV, and A con ^i 100 — 300 MeV, both 
very close to our Pade QCD findings. 

10. In their own novel gauge invariant nonperturbative approach to QCD 
at large distances, Gogohia and co-workers jP[| found that a momen- 
tum free parameter k in their approach, defined as some character- 
istic scale at which confinement and other nonperturbative {D\SB) 
effects begin to play a dominant role in QCD, or more broadly de- 
fined as a momentum that separates the nonperturbative phase of 
QCD from the perturbative phase, always has both a lower and an 
upper boundary value which Gogohia et. al. found numerically to be: 
635MeV < k < 775MeV. Their more recent estimates |45|, [46] give 
a higher value A x « 1.2 GeV, depending on input value for the pion 
decay constant. These result can be interpreted as a finding of a two 
scale structure in long distance QCD dynamics, and the explicit scale 
values estimated above by Gogohia et. al. are not incomparable with 
our paired scales shown in Tables |5| and ^| from Pade QCD. 

11. From a completely different phenomenological perspective based on 



known masses of hadronic resonances, Greiner and Schafer [47] had 
also argued that the transition from NPQCD to PQCD must occur 
rapidly within a small segment of the energy scale they found to be 
1 < VQ 1 < 3 GeV such that NPQCD is definitely all of Q < 1 GeV, 
while all of \[Cp > 3 GeV is PQCD. Stated in our own terms, this 
Greiner-Schafer argument places the onset point of NPQCD / DxS B 
in long distance QCD dynamics, in the range 1 < A x < 3 GeV, or 
1.0 < A x < 1.8 GeV, to be compared with the values: 0.706 < A x < 2.6 
GeV for Nf < 6, we found from our Pade QCD and shown variously in 
Tables [| and |]. 

12. On the specific issue of the confinement scale in QCD, Brown and 
Pennington [Q obtained an infrared plot based on considerations of 
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the Schwinger-Dyson equation for the gluon propagator, showing that 
the gluon propagator gets enormously enhanced and rises sharply in 
value to infinity in a Landau type cut-off behavior, at a scale Q ~ Aqcd 
in the infrared QCD region. This infrared gluon propagator behavior 
is normally taken as indicative of the point of onset of confinement in 
QCD. Its phenomenological scale found here at Q = Q c ps Aqcd is 
to be compared with our a^ or Y 2 plots shown variously in figs. §, |||, 
and figs. |6| to |9], which rise sharly to infinity at exactly similar critical 
momentum points Q C {Y 2 ) ~ 296MeV ~ Aqcd- 

13. The whole question of the enhanced Q~ A infrared singularity behavior 
of the gluon propagator has actually become a central pivot for formu- 
lating chiral symmetry breaking and confinement dynamics of QCD, 
and estimating the the scale boundaries of these phenomena. Gogohia 
and Magradze ESI used it to establish a close interrelatedness and jux- 



taposed existence in QCD of dynamical chiral symmetry breaking and 
confinement, exactly as our Pade QCD finds in figs. [I] to [^ and |6| to [l|. 
Arbuzov |){J working from the same premise of Q~ 4 infrared singularity 
behavior of the gluon propagator, derives various essential attributes 
of spontaneous chiral symmetry breaking and other nonperturbative 
infrared QCD phenomena, and finds a momentum scale k separating 
asymptotically free ultra violet QCD from the nonperturbative infrared 
domain to be k = 700 MeV for the specific flavor state Nf = he con- 
sidered. This Arbuzov value compares again very well with our value 
A x = 706 MeV for the same Nj = 0, in the R e+e -(Q) values shown in 
Table ||. Other recent studies |46], |5l], ^2| of the nonperturbative QCD 
vacuum, assumed dominated by a Q~ 4 gluon propagator behavior, are 
found to exhibit a two scale structure of chiral symmetry breaking and 
quark confinement dynamics, all strikingly similar to our Pade infrared 
QCD findings reported here. 



14. Finally Iwasaki and co-workers [Q found from their Lattice QCD stud- 
ies that dynamical chiral symmetry breaking and confinement in in- 
frared QCD must be closely related or correlated, even when they found 
no evidence that the two phenomena necessarily occur simultaneously. 
The implication for scale considerations, is that the two phenomena 
are compatible with being intrinsic two scale structure or sequential 
phenomena of infrared QCD dynamics, exactly as we find here with 
infrared Pade QCD. What is more, Iwasaki et. a!, found that their two 
scales and phenomena of infrared QCD, depend on the flavor number 
Nj of the QCD system, again exactly as our Tables ^ and ^|, and our 
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various Pade QCD plots, explicitly show. 

Clearly, these varied physical QCD phenomenological results concerning 
infrared QCD dynamics and its critical boundary points, all consistently close 
to and in agreement with the Pade QCD findings shown variously in figs, p] 
to H and figs. |6] to 10[ as well as Tables [5] and ||, compel us to admit the 



Pade infrared (03,04) structures as relating to physical QCD. In that case, 
we can make a number of important deductions from our work, concerning 
infrared dynamics of physical QCD as follows. 

6 FAVORED FLAVOR STATES FOR QUARK 
CONFINEMENT IN QCD 

One immediate deduction or proposal we can make following from the above, 
is that confinement and chiral symmetry breaking in QCD are caused by two 
distinct component color forces or dynamics of QCD. Of our two nonper- 
turbative infrared couplant forces Y 2 and Y"i, the latter is responsible for 
and identifiable with the onset and dynamics of chiral symmetry breaking 
in QCD, while the former Y 2 is responsible for quark confinement in QCD. 
Following from the fact that the two component color forces are indepen- 
dent solutions 03 and 04 of our Pade couplant equation, we assert as part 
of our finding, that the two nonperturbative infrared QCD phenomena, chi- 
ral symmetry breaking and quark confinement, are actually distinct physical 
phenomena. They may however still be closely related in the sense that one 
phenomenon provides a stimulus or background that facilitates the occur- 
rence of the other phenomenon. 

The indication from our plots is that dynamical chiral symmetry break- 
ing (DxSB) occurs first at the higher infrared nonperturbative scale Q c (Yi), 
phenomenologically flooding the QCD vacuum with qq pair condensates, 
which in turn provide a superconducting state of QCD vacuum, compati- 
ble with quark confinement in the manner advocated by a number of work- 
ers |54], |55], [56| . However, what our work now indicates is that confinement 
itself does not actually occur (even after a DxSB) until the independent Y 2 
confinement causing color force has attained its own confinement scale at 
Q = QdX'i) < QciXi)- in such a scenario, the possibility arises that if these 
two scales are very far apart, the condensate background effect caused by Y\ 
at scale Q c (Yi) may become dissipated or ineffective by the time we reach the 
Y 2 confining scale at Q C {Y 2 ). Then confinement does not occur even though 
we are at scale Q C (Y 2 ). 
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The details of this confinement mechanism do not concern us here. What 
we show now is that our work provides a strong evidence or indication that 
the above interrelatedness or correlation of confinement with dynamical chiral 
symmetry breaking does exist in concrete terms in QCD. Here we return to 
Table || and || where we already computed the ratios Q c (y 2 ) /Q c (Yi) as a 
measure of how close together or far apart the two nonperturbative infrared 
scales are. We now plot this Q c iX2)/QciXi) correlation factor as a function 
of flavor number and for different spacelike and timelike observables. The 
plots are shown in fig. |H| and reveal some interesting features. 

We observe that in the spacelike cases the correlation ratio first rises 
with increasing flavor number from Nf = 0, and reaches a peak at Nf = 2 
or in the region 1 < Nf < 3. It then falls off rapidly for all Nf > 4 and is 
practically zero near Nf = 7, 8. It is well known in nature, that the most 
stable hadrons and confined quark states are the (u, d) or Nf = 2 states 
of QCD. These are the nucleons (protons and neutrons). Nf = 3 confined 
quark states exist as higher baryon states that are only quasi stable. There 
are no known stable or quasi stable confined quark states of flavor number 
Nf > 4. Our fig. [13] appears therefore to be providing us a strikingly direct 



explanation and insight into the origin of this well known state of affairs 
in nature, attributing it to an interplay or some intrinsic interrelatedness of 
dynamical chiral symmetry breaking and quark confinement mechanism in 
QCD. 

In the case of the timelike Q cases shown also in fig. O, the correlation 



factor falls off progressively from Nf = to Nf = 8, in a manner to suggest 
that the most readily confined QCD states are the pure gluonic states or glue 
balls, rather than the 1 < Nf < 3 states. This would appear not to be so in 
nature but we are not able at the moment to reconcile this crucial difference 
in the timelike and spacelike behaviors of fig. [13|, except that the difference 
is directly traceable to p 2 being positive or negative in the manner shown in 
Table [TJ. Apart from this difference in their Nf < 3 behavior, the timelike and 
spacelike observables are fully in agreement from their correlation plots shown 
in fig. [IS], that confinement becomes rapidly less probable for Nf > 4, and 
that this decline arises from the intrinsic confinement scale Q c {Yi) becoming 
progressively less correlated with the intrinsic dynamical chiral symmetry 
breaking scale Q c (Yi) of QCD. 

We remark that Iwasaki and co-workers |J3|| actually found in their Lattice 



QCD studies, that confinement and dynamical chiral symmetry breaking 



cease to occur for Nf > 7, but occur only for Nf < 6. Our work and fig. [13 
explain adequately what Iwasaki and co-workers observed in their Lattice 
QCD studies, and this close agreement between their work and ours, leads 
us once more to affirm that our Pade infrared QCD results have a physical 
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reality, reproducing quite well infrared QCD dynamics from a wide variety 
of phenomenological considerations, including Lattice QCD simulations. 

7 FROZEN COUPLANT VERSUS INFRARED 
SLAVERY STRUCTURE OF CQCD 

A second deduction for QCD we can make from our Pade plots and features, 
is that while many of the phenomenological models we cited earlier in this 
paper, especially those dealing with ladder approximant Schwinger- Dyson 
equations, make mere guesses or assumptions about the form of the QCD 
coupling constant inside the CQCD regime, (A con < Q < A x ), some authors 
parameterizing it by a frozen coupling constant or a step function: 

9 2 ((q ~ P?) = 0(q 2 - p 2 )g 2 (q 2 ) + 9(p 2 - q 2 )g 2 (p 2 ) (35) 

with 

6(q 2 -p 2 ) = l,iiq 2 >p 2 

= 0, if q 2 <p 2 (36) 

or in the form: 
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where q c is some constant parameter identifiable with our A x , our Pade 
model shows explicitly [K| that the Nj < 8 flavor states are not infrared 
couplant freezing states of QCD. Rather the two coupling constant solutions 
a^ and a^ that govern QCD dynamics in this CQCD infrared regime of QCD 
are together growing infinitely large in the domain. Such infinitely growing 
coupling constant in the infrared region of QCD, is what is naturally needed 
to explain physical confinement of colored quarks and gluons into hadrons. 
Therefore, our optimized Pade QCD has this additional positive side to it 
that it offers a more natural explanation of the origin and mechanism of 
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quark confinement than some of the phenomenological models we cited in 
section 5. Indeed the Q~ A enhanced infrared singularity structure of the gluon 
propagator, currently gaining ground j|B], |5D], ^T], Q as a necessary ingredient 



to reconcile QCD vacuum with quark confinement, the linearly rising inter- 
quark potential law, and other nonperturbative features of infrared QCD, has 
no place for an infrared fixed point in low flavor states of QCD. The same 
conclusion that low flavored states of QCD do not manifest any infrared 
stable fixed point is found in several other places Refs. [[57]] to |66|], all rather 



contrary to Refs. |67|, |68 



8 SUMMARY AND CONCLUSIONS 

Summarizing, we state that the two scale cut off momentum structure exhib- 
ited in the infrared region of QCD by our optimized Pade couplant eqns. ([[]) 
and (§), is identifiable as marking the points of onset of dynamical chiral 
symmetry breaking and quark confinement respectively, in infrared QCD. 
Taking these two infrared nonperturbative phenomena of QCD as charac- 
teristically marking the beginning phase and end phase respectively, of the 
CQCD regime of QCD, we can conclude, based on extensive correlations 
of our results with a wide range of known phenomenological studies of low 
energy QCD, that our work actually constitutes a finding from within per- 
turbative QCD approaches, of a previously unknown but long speculated 
infrared QCD window, marked at its two ends by a critical coupling con- 
stant and a critical momentum, and separating QCD dynamics in the cold 
(zero temperature vacuum state), into the three distinct phases of purely 
hadronic QCD at very low momentum transfers Q, a largely weak perturba- 
tive QCD phase at high Q, and finally an intermediate phase or regime of 
strongly coupled yet quark-gluon dominated QCD dynamics. 

Traditionally in perturbative QCD (PQCD), one believes that an impene- 
trable wall at Q = Aq CD shuts out nonperturbative QCD and hadron physics 
from any view or access from the PQCD end. This apart, it is held that much 
of PQCD remains accessible in the domain stretch from the far ultraviolet 
(UV) QCD, until close to the Landau cut off wall at Q = Aqcd- What our 
work shows is that while the Landau wall at Q = Q C (Y2) ~ 296MeV ~ Aqcd 
remains impenetrable from the PQCD end, much of the most important as- 
pects of nonperturbative QCD dynamics including spontaneous chiral sym- 
metry breaking and the confinement of quarks into hadrons, actually lies in 
a buffer zone (CQCD) that leads up from the far UV QCD to the Landau 
wall. Beyond the Landau wall lies pure hadron physics, and our work pro- 
vides explicit boundary values both in critical momentum and critical QCD 
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coupling constant, separating these three phases of QCD, as delineated from 
the PQCD beta function parameters and renormalization group equations of 
a running QCD coupling constant. 

We believe the same CQCD buffer zone just delineated, is what is cur- 
rently being observed as Quark Gluon Plasma (QGP) from the hadron end 
at high temperature (non- vacuum QCD state), in heavy ion experiments 
at CERN and Brookhaven (RHIC). As such, it becomes immediately pos- 
sible to ask how the domain boundaries and dynamical characteristics of 
this intermediate CQCD/QGP phase of QCD get modified or shifted at two 
widely separated temperatures and energy density, represented by the heavy 
ion experiments on the one hand, and our zero temperature vacuum QCD 
studies on the other hand. This comparison will await further developments 
and evolution in the QGP plasma experiments and the various Lattice and 
phenomenological computations relating thereto. 

Focusing on our own boundary values of the CQCD regime computed 
here at zero temperature, we found that these values in general vary with 
flavor as seen from Tables |5| and ||, and from our various figures and plots. We 
have interpreted this Q c (Y2)/Qc(Xi) fl avor behavior to be a maifestation of 
a phenomenologically known fact that dynamical chiral symmetry breaking 
and quark confinement are closely related infrared QCD phenomena. Our 
specific finding is that the more widely separated their set scales are, the 
less likely the corresponding flavor state of QCD will be a quark confinement 
state. We actually found evidence from our fig. [IS] that quark confinement 
occurs mainly in the low flavor states Nf < 3, with a peak at Nf = 2, but 
becomes progressively non-existent or improbable in the higher flavor states 
Nf = 4, 5, 6, 7, 8. It turns out that this is actually what one finds in nature. 

After we had completed this work, our attention was drawn to a preprint 
by Elias and co-workers p9| in which they used their own method of de- 



nominator versus numerator zeros of Pade QCD beta function, to estimate 
the location and scale value of the point where perturbative QCD cuts off 
from the infrared region. They found only one scale cut-off point at a crit- 
ical momentum \x c « 1.0 GeV. This is in contrast to our two scale cut off 
structure presented variously in this paper, and anticipated in the several 
phenomenological models we cited in section 5 of this paper. We make two 
comments in relation to this work of Elias et. al. The first is that the method 
of denominator versus numerator zero of Pade beta function appears not ad- 
equate to reveal the multicomponent couplant solutions of Pade QCD and 
consequently the two scale infrared structure. This fact was already high- 
lighted in our earlier paper |TD[. One needs to work with effective charges 



and renormalization group invariants or optimized Pade QCD. The second 
comment we make is one of noting that the one scale cut off found by Elias 
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et. al. at fi c « 1.0 GeV, is comparable to our chiral symmetry breaking 
scale A x « 1.44 GeV. Since Elias et. al. worked with the higher Pade ap- 
proximants [2|1], [1|2], [3 1 1] , [1 1 3] , [2|2], as against our [1|1] Pade approximant, 
the close agreement between our onset scale and theirs, provides a good in- 
dication that our results and features presented here will hold irrespective 
of the particular Pade approximant used, exactly as we already inferred in 
paper fH|. One could of course repeat directly the various computations 
reported in this paper, using this time, these higher Pade approximants to 
verify that the results remain robust. The draw back is that the higher renor- 
malization scheme invariants p 3 ,p^... needed for such computations, are not 
presently known or computed perturbatively for the various QCD physical 
observables reported in this paper. 
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Table 1: Computed p2(Nf) values for various spacelike and timelike QCD 
observables 



Flavor 


Bj(e) 

H2 


Bj(u) 
P2 


GLS(u) 
H2 


(e+e") 
P2 


(t) 

P2 


(H) 

P2 


Number 
N f 

















17.9244 


15.9567 


17.9244 


-8.410032589173554 


1.8089 


-17.8458 


1 


13.7474 


12.3442 


13.3343 


-9.996607149709793 


-0.5654 


-18.9820 


2 


9.6051 


8.7025 


8.7787 


-10.91120013471925 


-2.9879 


-20.1684 


3 


5.4757 


5.0067 


4.2361 


-12.20710268197531 


-5.4756 


-21.4108 


4 


1.3304 


1.2238 


-0.3223 


-13.90995802777778 


-8.0508 


-22.7170 


5 


-2.8695 


-2.6906 


-4.9354 


-15.49181836878120 


-10.7438 


-24.0975 


6 


-7.1775 


-6.7981 


-9.6566 


-17.66469557734694 


-13.5963 


-25.5668 


7 


-11.6698 


-11.1857 


-14.5620 


-19.78668878025239 


-16.6675 


- 27.1456 


8 


-16.4581 


-15.9819 


-19.7635 


-22.74511792421761 


-20.0446 


-28.8640 


9 


-21.7139 


-21.3833 


-25.4325 


-25.96983971428571 


-23.8613 


-30.7687 


10 


-27.7143 


-27.7079 


-31.8461 


-30.64825148592373 


-28.3334 


-32.9345 


11 


-34.9379 


-35.5045 


-39.4829 


-36.70527878387512 


-33.8340 


-35.4912 


12 


-44.2885 


-45.8101 


-49.2467 


-46.58505774333333 


-41.0672 


- 38.6836 


13 


-57.7031 


-60.8467 


-63.0745 


-63.56012999671129 


-51.5427 


-43.0365 


14 


-80.2178 


-86.3914 


-86.0023 


-101.9145678055556 


-69.1847 


-49.9139 


15 


-130.2979 


-143.6271 


-136.4956 


-229.8874551066667 


-108.3829 


-64.3179 


16 


-374.1328 


-423.2004 


-380.7437 


-1724.404563921111 


-298.6440 


-131.2671 
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Computed P2JNf) values for various spacelike and timelike 



QCD observables; 33 



2 Yi structural parameters Q m j n , ai(Qmm), Qmax and a^jQi 



for Bjorken sum rule /y values shown in Table 1| in opti- 
mized [1|1] Pade QCD. The Nf = 0, 1, 2 cases were computed 



using 
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3 Bjorken neutrino sum rule Y"x structural parameters Q m j n , ai(Qmm), Q 



and a3{Q max ) computed using the p? values shown in Ta- 



WTj, in optimized [1|1] Pade QCD. The N f 
were computed using J\D\ 



0, 1, 2 cases 
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4 Tau(r) Hadronic decay Y 1 structural parameters Q m j n , aiiQmm), Q 



and Q3(Qmax) computed using the pi values shown in Table 1 
in optimized [1|1] Pade QCD. All D values in this case, in- 
cluding Nf = state, are positive 
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Table 2: Yi structural parameters Q min , ai(Q min ), Q max and a 3 (Q ma , x ) for 



Bjorken sum rule p 2 values shown in Table [TJ in optimized [1|1] Pade 
QCD. The Nf = 0,1,2 cases were computed using 



Flavor Number Nf 


Qmin (GeV) 


O'liQmm) 


<5ma X (GeV) 


a 3(vmax) 





16.3391 


0.0812 


633.07 


0.10006 


1 


12.5361 


0.0936 


444.10 


0.1174 


2 


8.0622 


0.1164 


325.33 


0.15139 


3 


6.11284 


0.13944 


364.81 


0.18538 


4 


11.1934 


0.12027 


238.4 


0.1477 


5 


17.012 


0.1149 


204.84 


0.1344 


6 


24.2354 


0.11425 


161.30 


0.1278 


7 


32.592 


0.11665 


115.8 


0.1243 


8 


36.214 


0.12191 


37.56 


0.12235 


9 


6.910 


0.12087 


36.279 


0.1303 


10 


0.3015 


0.11868 


62.46 


0.1428 


11 


1.036(10~ 3 ) 


0.1143 


250.952 


0.16116 


12 


8.2160(10~ 7 ) 


0.1059 


7.7850(10 +3 ) 


0.1900 


13 


9.029(10" 14 ) 


0.0918 


7.8324(10 +7 ) 


0.2395 


14 


9.752(10~ 35 ) 


0.07130 


7.8451(10+ 20 ) 


0.3652 


15 


1.281(KT 101 ) 


0.0448 


w 1.0(10+ 85 ) 


>0.80 


16 


2.250(10" 176 ) 


0.01520 


> 10+ 308 


> 1.0 
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sum rule Y\ structural parameters 



Table 3: Bjorken neutrino 

Vmin) a l(vmin)) Vmax an d a 3\Q\ 

shown in Table [l], in optimized [1|1] Pade QCD. The Nf 
computed using \I\D 



computed using the p 2 values 

0,1,2 cases were 



Flavor Number Nf 


Qmin (GeV) 


^1 (.^imin) 


Qmax (GeV) 


a 3(vmax) 





16.271 


0.0868 


404.4 


0.1087 


1 


12.432 


0.1005 


403.6 


0.1284 


2 


7.63 


0.1285 


326.6 


0.1724 


3 


8.5585 


0.1312 


266.8 


0.171 


4 


13.5895 


0.1194 


732.4 


0.1464 


5 


19.388 


0.1159 


281.4 


0.1357 


6 


26.4586 


0.1159 


343.3 


0.1299 


7 


34.4886 


0.1186 


113.8 


0.1265 


8 


37.52 


0.1237 


39.72 


0.12415 


9 


5.117 


0.1219 


37.18 


0.1315 


10 


0.2723 


0.1187 


64.03 


0.1430 


11 


1.678(10- 3 ) 


0.1133 


260.28 


0.1590 


12 


8.130(10- 8 ) 


0.10402 


8.2617(10+ 3 ) 


0.1835 


13 


9.620(10~ 15 ) 


0.0897 


8.592(10+ 7 ) 


0.2255 


14 


6.69(10- 42 ) 


0.06957 


8.9984(10+ 20 ) 


0.3253 


15 


3.440(10~ 83 ) 


0.0441 


1.02687(10+ 85 ) 


1.465 


16 


4.90(10- 307 ) 


0.015040 


> 10+ 305 


>10.0 
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Table 4: Tau(r) Hadronic decay 

<5min,oi(<3min),Qmax and a 3 (Q max ) computed 
Table 0, in optimized [1|1] Pade QCD. All D 
Nf = state, are positive. 



Y\ structural parameters 
using the p\ values shown in 
values in this case, including 



Flavor Number Nf 


Qmin (GeV) 


ai(Qmin) 


Qmax (GeV) 


6^3 1^^; max J 





26.6430 


0.0773 


765.30 


0.0942 


1 


32.0952 


0.07825 


1697.95 


0.0942 


2 


38.7402 


0.07977 


749.06 


0.09475 


3 


46.851 


0.08194 


785.50 


0.0959 


4 


56.716 


0.0849 


813.4 


0.09772 


5 


68.5620 


0.08893 


660.9 


0.10015 


6 


82.1620 


0.09422 


547.18 


0.1033 


7 


95.261 


0.10123 


267.0 


0.10692 


8 


92.2863 


0.11058 


95.54 


0.11094 


9 


15.2460 


0.1147 


81.01 


0.1232 


10 


0.591 


0.1172 


122.64 


0.1407 


11 


3.5330(10~ 3 ) 


0.1164 


435.71 


0.1654 


12 


1.2110(10- 6 ) 


0.11016 


1.2066(10+ 4 ) 


0.2035 


13 


2.015(10- 14 ) 


0.0964 


1.09640(10+ 8 ) 


0.2728 


14 


3.505(10- 33 ) 


0.0745 


1.00163(10+ 21 ) 


0.4690 


15 


2.636(10~ 105 ) 


0.0463 


> 1.0(10+ 88 ) 


>1.0 


16 


1.0610- 295 


0.0153 


> 10+ 305 


>1.0 
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Table 5: CQCD boundary momenta Q C (Y\_) and Q C (Y 2 ) computed for var- 
ious physical observables or effective charges in the [1|1] optimized Pade 
QCD. The Q C {Y X ), Q C {Y 2 ) values are taken from figs, g |, and figs. § to 0, 
at the three strategic couplant points discussed in the paper. Ratio I = 
Qon(Y 2 )/Qon(Yi), while Ratio III = Q ff(Y 2 )/Qoff(Yi). Ratio II was evalu- 
ated at the IR attractor point of the Y 2 curve where this exists. 



QCD 
observable 


Flavor 
N, 


Qon{Y 2 ) 

GeV 


Qon\Y\) 

GeV 


Ratio 
I 


Qoff(Y 2 ) 
GeV 


Qoff(Yi) 

GeV 


Ratio 
III/II 


Bj(e) 





1.0 


2.375 


0.4211 


0.8791 


0.9465 


0.9287 




1 


1.195 


2.50 


0.4780 


1.0956 


0.900 


1.2173 




2 


1.60 


2.825 


0.5664 


1.5661 


0.8380 


1.8689 


! 


3 


1.815 


3.25 


0.5585 


0.85 


0.81 


1.0494 




4 


1.09 


3.186 


0.3421 


0.78278 


0.88625 


0.8833 




5 


0.8125 


3.60 


0.2257 


0.3953 


0.98629 


0.4008 




6 


0.6340 


4.25 


0.1492 


0.2450 


1.1392 


0.2150 




8 


0.440 


7.80 


0.0564 


0.1084 


1.900 


0.0571 


BjM 





1.4474 


2.9460 


0.4913 


1.349312 


1.9650 


0.6867 




1 


1.7177 


3.10 


0.5541 


1.66233 


2.3360 


0.7116 




2 


2.4460 


3.8150 


0.6412 


- 


- 


- 




3 


1.9365 


3.760 


0.5150 


1.856 


2.5780 


0.7199 




4 


1.30 


3.910 


0.3325 


0.911941 


1.3602 


0.6704 




5 


0.95 


4.10 


0.2317 


0.4732 


1.15174 


0.4109 




6 


0.75 


4.870 


0.1540 


0.290 


1.2941 


0.2241 




8 


0.4710 


8.0 


0.0589 


0.1225 


2.0593 


0.0595 


Re+e-{Q) 




1 


0.3700 
0.350 


1.700 

1.850 


0.2176 
0.1892 


0.2364 
0.20996 


0.7064 
0.7345 


0.3347 

0.2859 




2 


0.337 


2.05 


0.1644 


0.19005 


0.7611 


0.2497 




3 


0.320 


2.175 


0.1471 


0.16717 


0.80043 


0.2089 




4 


0.3037 


2.625 


0.1157 


0.1425 


0.8599 


0.1657 




5 


0.2875 


3.00 


0.0958 


0.1204 


0.9420 


0.1278 




6 


0.270 


3.800 


0.0711 


0.09771 


1.07991 


0.0905 




7 


0.250 


4.80 


0.0521 


0.079425 


1.320 


0.0602 




8 


0.249 


7.25 


0.0343 


0.06250 


1.8115 


0.0345 


t decay 





1.2446 


3.373 


0.369 


1.07444 


1.340 


0.8018 


R& 


1 


1.1341 


3.674 


0.3087 


0.8308 


1.3980 


0.5943 




2 


1.025 


4.0 


0.2562 


0.6615 


1.468 


0.4506 




3 


0.9340 


4.50 


0.2076 


0.5305 


1.560 


0.3401 




4 


0.850 


5.237 


0.1623 


0.42325 


1.6825 


0.2516 




5 


0.7756 


6.250 


0.1241 


0.3331 


1.8540 


0.1796 




6 


0.700 


7.50 


0.0933 


0.260 


2.146 


0.1212 




8 


0.6170 


14.70 
1 Ll 37— 


0.0420 


0.15226 


3.590 


0.0424 



Table 6: Further data on CQCD boundary momenta Q c (Yi) and Q C (Y 2 ) 
computed for various physical observables or effective charges in the [1|1] 
optimized Pade QCD. The Q c (Yi),Q c (Y 2 ) values are taken from figs. ||, |3|, 
and figs. ^] to |T2|, at the three strategic couplant points discussed in the paper. 
Ratio I = Qon(Y 2 )/Q n(Y 1 ), while Ratio III = Q ff(Y 2 )/Q off (Y 1 ). Ratio II 
was evaluated at the IR attractor point of the Y 2 curve where this exists. 
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Y1 Profile plots 
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ure 4: A plot showing similar Y\ triple point crossing feature or solutions 
a2, a?, for various physical observables, spacelike and timelike, exactly as 
the R e +e- (Q) observable discussed in paper fTTJ . 
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Y2 Profile plots 
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Figure 5: A plot showing similar one crossing point feature or solution 0,4 of 
Y 2 for various QCD observables, spacelike and timelike, except for Nf < 3 in 
the spacelike cases, when Y" 2 exhibits a double or triple crossing point feature. 
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Figure 6: A two scale cut off structure of 03 and 04 in the infrared region: 
The case of Bjorken polarized sum rule Bj(e), (spacelike). 
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Figure 7: A two scale cut off structure of 03 and 04 in the infrared region: 
The case of Bjorken neutrino sum rule Bj(z/), (spacelike). 
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IR cut off points of a3,a4 with Tau hadronic decay coefficients 
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Figure 8: A two scale cut off structure of 03 and GI4 in the infrared region: 
The case of tau hadronic decay rate, (timelike). 
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IR cut off points of a3,a4 with Higgs decay width coefficients 
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Figure 9: A two scale cut off structure of a^ and 04 in the infrared region: 
The case of Higgs hadronic decay width, (timelike). 
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IR cutoffs of Y1 andY2, Nf=2 
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Figure 10: The Nf = 2 infrared cut offs of Y-± and Y 2 plotted for two spacelike 
observables for which p 2 is positive and large, and D is negative. The cut off 
curves are seen to bunch closely together. 
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Figure 11: The cases of Nf — 0, 1, infrared cut 
the two spacelike observables, Bj(e) and Bj(z/) 
large, and D is negative 
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Figure 12: A plot showing the Y 2 cut off behavior in the Nf = 3,4 flavor 
states of positive p 2 but positive D, in spacelike physical observables, Bj(e), 
Bj(z/), and Gross - Llewellyn Smith (GLS) sum rules. The cut offs a^,a^ are 
seen to occur closely together, with a^ even growing faster than a% in the 
ensuing CQCD dynamics in Nf = 3 case 
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Y1 ,Y2 correlation in IR Pade QCD 
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Figure 13: A plot showing the strong correlation between dynamical chiral 
symmetry breaking and quark confinement in QCD as revealed in Pade in- 
frared couplant structures. Quark confinement is seen to be most favored in 
Nf < 3 states with a peak at Nf = 2, but drops rapidly to zero for Nf > 4, 
exactly as one finds in nature. Each curve was normalized to its Nf = 
flavor state. 
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5 CQCD boundary momenta Q C {Y\) and Q C {Y 2 ) computed for 



various physical observables or effective charges in the 1 1 1 1 1 op- 
timized Fade QCD. The Q c {Yi), Qc{Y 2 ) values are taken from| 



figs. 2| , ||, and figs. |6] to [12], at the three strategic couplant 
points discussed in the paper. Ratio I = Q an ,(Y2)/Qo n (Yi), 
while Ratio III = Qoff(Y2)/Qoff(Yi). Ratio II was evaluated 
at the IR attractor point of the Y 2 curve where this exists. . 



6 Further data on CQCD boundary momenta Q C {Y\) and Q c ( Y 2 
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computed for various physical observables or effective charges 



in the [1|1J optimized Fade C)CD. The Q c {Yi),Q c {Y 2 ) val- 



ues are taken from figs.~2| , |3], and figs. ^| to 12, at the three 
strategic couplant points discussed in the paper. Ratio I = 
Qon(Y 2 )/Q 0n (Y 1 ), while Ratio III = Q off (Y 2 )/Q off (Y 1 ). Ra- 
tio II was evaluated at the IR attractor point of the Y 2 curve 
where this exists 
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